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ABSTRACT 

Thomson optical depth r measurements from Planck provide new insights into the reionization of 
the universe. In pursuit of mo del-independent constraints on the properties of the ionizing sources, we 
determine the empirical evolution of the cosmic ionizing emissivity. We use a simple two-parameter 
model to map out the evolution in the emissivity at z >6 from the new Planck optical depth r 
measurements, from the constraints provided by quasar absorption spectra and from the prevalence of 
Lyman a emission in z ^ 7-8 galaxies. We find the redshift evolution in the emissivity Nion{z) required 
by the observations to be (dlog^g Von/d^)2=8 = —0.15lo °® ((d log^g iVion/d2:)2=8 = —0.19toui for a 
flat prior), largely independent of the assumed clumping factor Chii and entirely independent of the 
nature of the ionizing sources. The trend in Nio^{z) is well-matched by the evolution of the galaxy 
C/U-luminosity density (dlog^g Puv /dz = —0.11 ± 0.04) to a magnitude limit > —13 mag, suggesting 
that galaxies are the sources that drive the reionization of the universe. The role of galaxies is 
further strengthened by the conversion from the UV luminosity density puv to Nion{z) being possible 
for physically-plausible values of the escape fraction fesc^ the Lyman-continuum photon production 
efficiency ^lon, and faint-end cut-off Mum to the luminosity function. Quasars/AGN appear to match 
neither the redshift evolution nor normalization of the ionizing emissivity. Based on the inferred 
evolution in the ionizing emissivity, we estimate that the z ~ 10 UV-luminosity density is 8/^4® x 
lower than at z ^ 6, consistent with the observations. The present approach of contrasting the 
inferred evolution of the ionizing emissivity with that of the galaxy UV luminosity density adds to the 
growing observational evidence that faint, star-forming galaxies drive the reionization of the universe. 
Subject headings: galaxies: evolution — galaxies: high-redshift 


1. INTRODUCTION 

One of the most important phase transitions in the 
history of the universe is the reionization of the neu¬ 
tral hydrogen gas. Following recombination early in 
the universe shortly after the Big Bang, the universe 
likely remained in a largely neutral state until z ^ 15- 
25. The collapse of the first dark matter halos and gas 
cooling brought about the formation of the first stars 
and galaxies. These early stars and galaxies have long 
been thought to provide the ionizing photons necessary 
to reionize the universe (Loeb & Barkana 2001; Loeb 
2006). 

Despite this general picture of reionization and the 
likely role that early galaxies played in the process, es¬ 
tablishing that this is the case has been particularly chal¬ 
lenging, both due to the difficulties in probing the ioniza- 

^ Based on observations made with the NASA/ESA Hubble 
Space Telescope, which is operated by the Association of Uni¬ 
versities for Research in Astronomy, Inc., under NASA contract 
NAS 5-26555. 

^ Leiden Observatory, Leiden University, NL-2300 RA Leiden, 
Netherlands 

® UCO/Lick Observatory, University of California, Santa 
Cruz, CA 95064 

^ Department of Astronomy, Yale University, New Haven, CT 
06520 

® Leibniz-Institut fur Astrophysik Potsdam (AIP), An der 
Sternwarte 16, D-14482 Potsdam, Germany 

® Department of Physics and Astronomy, South Road, 
Durham, DHl SEE, United Kingdom 

^ Department of Physics & Astronomy, University of Sussex, 
Palmer, BRIGHTON, BNl 9QH, United Kingdom 


tion state of the universe at z > 6 (e.g., Ouchi et al. 2010; 
Stark et al. 2010) and providing constraints on the ion¬ 
izing photons that early galaxies themselves are thought 
to produce (e.g., Siana et al. 2010, 2015; Vanzella et al. 
2012; Nestor et al. 2013; Mostardi et al. 2013; Cooke et 
al. 2014). Further complicating the interpretation were 
early measurements of the Thomson optical depth by 
WMAP (e.g., T = 0.17 ±0.06: Spergel et al. 2003) which 
suggested a significant quantity of the ionizing photons in 
the early universe - implying an instantaneous reioniza¬ 
tion redshift of Zr = 20itg^ (95% confidence) - and which 
showed no clear connection with the early evolution of 
galaxies. While lower optical depths were measured in 
subsequent studies by WMAP (i.e., r = 0.089 ± 0.014: 
Bennett et al. 2013) with an implied z^ = 10.6 ± 1.1, 
these depths still pointed towards substantial amounts 
of ionizing radiation being present in the early universe. 

Fortunately, substantial progress has been made over 
the last ten years to better understand cosmic reion¬ 
ization. Much of the progress has been observational, 
through the better study of bright quasars and improved 
statistics on Lya emission in normal star-forming galax¬ 
ies, to better probe the ionization state of the z = 6- 
9 universe (e.g., Mortlock et al. 2011; Ono et al. 2012; 
Pentericci et al. 2014 [P14]; Schenker et al. 2014 [S14]; 
Caruana et al. 2014). The greater depths of probes for 
z = 6-10 galaxies provided greater confidence that galax¬ 
ies could provide the necessary reservoir of photons to 
reionize the universe (e.g., Oesch et al. 2010; Bunker et 
al. 2010; Bouwens et al. 2010, 2011; Ellis et al. 2013; 
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TABLE 1 

Key Observational Constraints on the Reionization History of the Universe.® 


Constraint 





# 

Redshift 

Qhii Constraint 

Technique 

Reference 


Key Constraints on the Ionization History of the Universe Explicitly Considered Here 



In Deriving the Evolution of the Cosmic Ionizing Emissivity at 2 > 6 

1. 


Thomson Optical Depth r = 0.066 zb 0.013 

PC15 

2. 


Reionization Finishes Between 2 : = 5.9 and 2 : = 6.5 



5.03 

u.y»yy‘±oj._Q QQQQjgg 

Gunn-Peterson Optical Depth 

Fan et al. (2006a) 


5.25 

u.yyy»oou_Q QQQQ244 

Gunn-Peterson Optical Depth 

Fan et al. (2006a) 


5.45 

0 9999333+0-0000247 
u.y»yyooo_Q QQQQ3Qj^ 

Gunn-Peterson Optical Depth 

Fan et al. (2006a) 


5.65 

QQQQ4gQ 

Gunn-Peterson Optical Depth 

Fan et al. (2006a) 


5.85 

0.99988001°°“ 

Gunn-Peterson Optical Depth 

Fan et al. (2006a) 


6.10 

0.99957±0.00030 

Gunn-Peterson Optical Depth 

Fan et al. (2006a) 


5.9 

>0.89 

Dark Gaps in Quasar Spectra 

McGreer et al. (2015) 


5.6 

>0.91 

Dark Gaps in Quasar Spectra 

McGreer et al. (2015) 


6.24-6.42 

<0.9 (2(t) 

Lyo: Damping Wing of Quasars 

Schroeder et al. (2013) 



Higher-Redshift Constraints 


3. 

7.0 

Qniiiz = 7) = omtoll 

Prevalence of Lyo Emission in Galaxies 

S14 

4. 

8.0 

Qhii{z = 8) < 0.35 

Prevalence of Lyo Emission in Galaxies 

S14 



Continuity with Ionizing Emissivity Estimates at 2 = 4.75 


5. 


logiQ Nion{z = 4.75) = 

10°°-99±0-45 s-1 Mpc-3 

BB13 


Other Constraints on the Ionization History of the Universe Not Explicitly Used^ 


6.3 

>0.5 

Lyo Damping Wing of a GRB 

Totani et al. (2006) 
McQuinn et al. (2008) 


6.6 

>0.6 

Lyo Emitters 

Ouchi et al. (2010) 


6.6 

>0.5 

Galaxy Clustering 

McQuinn et al. (2007), 





Ouchi et al. (2010) 


7.0 

0.32-0.64 

Lycf-Emitter LFs 

Ota et al. (2008) 


7.0 

~0.5 

Prevalence of Lyo Emission in Galaxies 

Caruana et al. (2014) 


7.0 

0.1-0.4 

Prevalence of Lyo Emission in Galaxies 

Ono et al. (2012) 


7.0 

<0.49 

Prevalence of Lyo Emission in Galaxies 

P14 


7.0 

<0.5 

Prevalence of Lyo Emission in Galaxies 

R13® 


7.0 

<0.5 

Clustering of Lyo Emitting Galaxies 

Sobacchi & Mesinger (2015) 


7.1 

<0.9 

Near-Zone Quasar 

Mortlock et al. (2011), 





Bolton et al. (2011) 


8.0 

<0.70 

Prevalence ot Lyo Emission in Galaxies 

Tilvi et al. (2014) 


® This table is a compilation of the constraints presented in the original papers under References, but with valuable guidance 
by the results presented in Figures 5 and 3 from R13 and R15, respectively. 

^ While not explicitly considered in deriving the evolution of the cosmic ionizing emissivity, almost all of these constraints are 
satisfied for the typical reionization histories derived in this study (see right panel of Figure [9t. 

R13 estimate this constraint on Qhii based on the observational results from Fontana et al. (20f0), Pentericci et ah (2011), 
Schenker et al. (2012), and Ono et al. (2012) and the simulation results from McQuinn et al. (2007), Mesinger h Furlanetto 
(2008), and Dijkstra et al. (2011). 


McLure et al. 2013; Oesch et al. 2013; Bouwens et al. 
2015 [B15]). However, theoretical progress has been sim¬ 
ilarly substantial, due to significantly improved estimates 
of the clumping factor (Bolton & Haehnelt 2007; Paw- 
lik et al. 2009, 2015; Finlator et al. 2012; Shull et al. 
2012 ) and ever more sophisticated simulations tracking 
the reionization of the universe and the propagation of 
Lya out of the galaxies and into the IGM (e.g., Mesinger 
et al. 2015; Choudhury et al. 2015). Lastly, measure¬ 
ments of the integrated column of ionized material to the 
last-scattering surface from 3-year Planck mission yield 
T — 0.066±0.013 (Planck Collaboration et al. 2015, XIII 
[PCI5]), implying that Zr = ^-nd suggesting that 

current surveys may be uncovering the sources that led 
to the reionization of the universe. 

Naturally, numerous studies have taken advantage of 
this collective progress to construct self-consistent mod¬ 
els for reionizing the universe (e.g., Choudhury & Fer¬ 


rara 2006; Bolton & Haehnelt 2007; Oesch et al. 2009; 
Trent! et al. 2010; Haardt & Madau et al. 2012 [HM12]; 
Bouwens et al. 2012a; Kuhlen & Faucher-Ciguere 2012 
[KF12]; Shull et al. 2012; Finkelstein et al. 2012b; Al¬ 
varez et al. 2012; Robertson et al. 2013 [R13]; Cai et al. 
2014; Choudhury et al. 2015; Ishigaki et al. 2015; Robert¬ 
son et al. 2015 [R15]), and it is indeed encouraging that 
many recent models (e.g., R15) prove successful in reion¬ 
izing the universe, while matching many other individual 
constraints on the reionization state of the universe and 
also the consensus star formation history (e.g., Madau & 
Dickinson 2014; B15). 

While these analyses are reassuring and offer strong 
support for the idea that galaxies reionize the universe, 
the uniqueness of galaxies as the source of photons to 
reionize the universe is more challenging to establish. 
Indeed, it is possible to imagine the existence of other 
populations of ionizing sources (e.g., Chen et al. 2003; 






3 



Redshift 

Fig. 1.— The key observational constraints on the filling factor 
of ionized hydrogen Qhii considered here in modeling the evo¬ 
lution of the cosmic ionizing emissivity (§2 and Table [Q. These 
constraints include the Gunn-Peterson optical depths and dark- 
gap statistics measured in 2 : ~ 6 quasars (Fan et al. 2006a; 
McGreer et al. 2015), damping wings measured in 2 : ~ 6.2-6.4 
quasars (Schroeder et al. 2013), and the prevalence of Lyo; emis¬ 
sion in 2 : ~ 7-8 galaxies (S14). Based on the constraints shown 
in this figure (particularly on the basis of the Gunn-Peterson 
troughs and dark-gap statistics and damping wings of 2 : ~ 6-6.5 
quasars), reionization ends between 2 = 5.9 and 2 = 6.5. The 
other key constraints we consider are the new Thomson optical 
depth measurements from Planck (PG15: represented here with a 
upward-pointing red arrow at the instantaneous reionization red- 
shift Zreion = ^ requirement for continuity with the 

cosmic ionizing emissivity at 2 = 4.75, as derived by BB13. The 
large changes in the filling factor are represented schematically by 
the grey-shaded region. 

Somerville et al. 2003; Hansen & Haiman 2004; Madau 
et al. 2004; Ricotti & Ostriker 2004; Ricotti et al. 2008), 
which even if speculative could also match current con¬ 
straints. 

What model-independent statements can be made 
about the sources that reionize the universe? To an¬ 
swer this question, it is useful to look at the evolution 
of the cosmic ionizing emissivity, since this allows us to 
keep assumptions regarding the nature of the ionizers to 
a minimum. As we will show, interesting constraints on 
the evolution of the emissivity can be obtained based 
on current observations, e.g., from the Thomson optical 
depths or the inferred filling factor of ionized hydrogen 
Qhii at z = 6-9 (see also Mitra et al. 2011, 2012, 2013). 

The purpose of the present analysis is to take advan¬ 
tage of current observational constraints on the ioniza¬ 
tion state of the universe at z > 6 (e.g.. Fan et al. 2006a; 
S14; PC15) to constrain the evolution of the cosmic ioniz¬ 
ing emissivity with redshift. Through comparisons with 
the evolution of the 17V luminosity density of galaxies 
and other potential ionizing sources, we can evaluate the 
likelihood that each of these sources of ionizing photons 
drives the reionization of the universe. We begin with a 
description of the relevant observations (§2) and method¬ 
ology (§3) and then derive constraints on the evolution of 
the ionizing emissivity Vion (§4.1). After deriving con¬ 
straints on Nio^(z)^ we compare our results with what 
we would expect for galaxies (§4.2), quasars (§4.3), and 


consider the implications for the UV luminosity density 
at 2 ~ 10 (§4.4). We conclude with a brief summary 
(§5). We take Hq = 67.51 ± 0.64, Ha = 0.6879 ± 0.0087, 
and = 0.3121 ± 0.0087, = 0.02230 ± 0.00014 

(PC15). 

2. OBSERVATIONAL CONSTRAINTS 

There are a wide variety of observational constraints in 
the literature on the ionization state of neutral hydrogen 
that can be leveraged in considering questions regarding 
the reionization of the universe. 

In the present analysis, we consider constraints at four 
different points in the reionization history of the universe. 
First of all, there are useful constraints on the end of 
cosmic reionization at z = 5.0-6.5. The most important 
constraints make use of the Gunn-Peterson optical depth 
measured from bright z & quasars (Fan et al. 2006a,b; 
Becker et al. 2001) and also by looking at the distribu¬ 
tion of dark gaps in x = 6 quasar spectra (McGreer et 
al. 2015) which suggest that cosmic reionization is com¬ 
plete by z = 5.9. By contrast, the presence of damp¬ 
ing absorption wings in the spectra of three z = 6.2-6.4 
quasars studied by Schroeder et al. (2013) strongly sug¬ 
gest that cosmic reionization is not complete by z = 6.4. 
These results clearly indicate that cosmic reionization is 
completed between z = 6.5 and z = 5.9. 

Second, we consider the constraints on the reionization 
history of the universe that come from the prevalence of 
Lya emission in the UV spectra of z ~ 6-8 galaxies. 
Assuming that the prevalence of Lya emission in star¬ 
forming galaxies at z 6.5-8 is a simple extrapolation of 
the observed prevalence at z ~ 4-6 and any departures 
from these trends are due to an increasingly neutral IGM 
at z ~ 7-8, one can use this technique to quantify the 
filling factor of ionized hydrogen at z ^ 6.5-8 (Santos 
2004; Malhotra & Rhoads 2004; McQuinn et al. 2007; 
Mesinger & Furlanetto 2008; Stark et al. 2010, 2011; 
Fontana et al. 2010; Dijkstra et al. 2011; Pentericci et 
al. 2011; Ono et al. 2012; Treu et al. 2013; Caruana 
et al. 2012, 2014; Tilvi et al. 2014; but see also Bolton 
& Baehnelt 2013; Mesinger et al. 2015). Here, we will 
make use of the recent constraints on Qhii from S14, 
i.e., Qhii{z = 7) = *3.nd Qhii{z = 8) < 0.35, 

using the McQuinn et al. (2007) models, but also briefly 
consider the impact of an alternate set of constraints. 

Of course, the real uncertainties on Qhii are likely 
larger than the formal uncertainties quoted by S14, as 
can be seen by comparing the different estimates for 
Qhii quoted by S14 based on the same set of observa¬ 
tions or using the results from other studies (Pentericci 
et al. 2014; Tilvi et al. 2014). Therefore, we take the 
formal Icr uncertainty in Qhii to be ±0.15, such that 
Qhii{z = 7) = 0.66 ± 0.15. Consistent with the con¬ 
straints from S14, we suppose that Qhii{z = 8) is equal 
to 0.15±0.20. However, we recognize that the constraints 
we use are just estimates and the true values could be 
different if the assumptions used in deriving these frac¬ 
tions are not correct (e.g., regarding the velocity offset 
for Lya line or the photoionization rate r(z): Mesinger 
et al. 2015). 

Third, we consider the constraints on the cosmic ioniz¬ 
ing emissivity that come from the Thomson optical depth 
r = 0.066 ± 0.013 measured from the three-year Planck 
results (PC15). While this optical depth measurement 
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Fig. 2.— Allowed evolution (68% confidence) of the cosmic ionizing emissivity over the redshift range z = 6 to z = 11, as inferred from 
the observations, {upper left) Constraints on the evolution of the emissivity supposing a completion of reionization at z = 5.9-6.5 (Fan 
et al. 2006a; Schroeder et al. 2013; McGreer et al. 2015) and the 9-year measured Thomson optical depth r = 0.089 ± 0.14 from WMAP 
(Bennett et al. 2013), and a simple modeling of the evolution of Qhii with Chii = 3. The region so defined is presented in hatched 
magenta. Constraints on the emissivity from z = 2 to z = 5 from BB13 are shown with the light-blue-shaded region and were derived from 
observations of the Lya forest. Systematic uncertainties in the inferred emissivity can often be large, ~ 0.3 dex, due to the sensitivity to 
the assumed (or measured) temperature of the IGM and the opacity to ionizing photons. As a result, some past constraints have been much 
lower than those from BB13 (e.g., KF12: dotted-grey-hatched region). The light gray line shows the ionizing emissivity model from HM12 
constructed to match the WMAP r’s (and the models from KF12 are similar to this [see the lower panel of Figure 10 from KF12]). {upper 
right) Identical to the upper-left panel except including the optical depth constraints from Planck (PC15: red-shaded region), {lower left) 
Constraints on the evolution of the ionizing emissivity {orange-shaded and green-shaded regions) assuming a completion of reionization at 
z = 5.9-6.5, the Planck Thomson optical depths, and various constraints on Qhii using the prevalence of Lyo emission in candidate z ~ 7-8 
galaxies, {lower right) Similar to the lower left panel, but using the Qhii constraints from S14 {red-shaded region) and also showing the 
95% confidence intervals. The 68% and 95% confidence intervals featured in this panel are also expressed in tabular form in Table [2] 
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TABLE 2 

68% AND 95% Confidence Intervals on 
THE Inferred Ionizing Emissivity 
VERSUS ReDSHIFT FOR A FIDUCIAL 

clumping factor ChII of 3.^ 


logiQ iVion [s ^ Mpc 3] 


Redshift 

Lower Bound 
95% 68 % 

Upper Bound 
68 % 95% 

6 

50.72 

50.81 

51.04 

51.11 

7 

50.69 

50.73 

50.85 

50.90 

8 

50.52 

50.60 

50.71 

50.74 

9 

50.27 

50.41 

50.62 

50.69 

10 

50.01 

50.21 

50.56 

50.66 

lO 

49.75 

50.00 

50.49 

50.64 

12 ^’ 

49.51 

49.80 

50.43 

50.61 

Id*’ 

49.24 

49.60 

50.36 

50.59 

14 b 

48.99 

49.39 

50.29 

50.57 

15*> 

48.74 

49.18 

50.23 

50.55 


® The results tabulated here are featured in 
the lower right panel of Figure [5] and make 
use of all 5 key observational constraints 
considered here (Table [!}• These results are 
derived in the context of the simple two- 
parameter model described in §3. 

^ Results here more sensitive to r’s mea¬ 
surements from Planck and functional form 
adopted in modeling the ionizing emissivity 
evolution. 

does not provide any information on the filling factor 
of ionized hydrogen at a specific time, it does provide 
a powerful constraint on the integrated path length of 
ionized hydrogen to the last-scattering surface. 

Finally, we require that the ionizing emissivity extrap¬ 
olated to z = 4.75 be consistent with the io50-99±0-45 g-i 
Mpc“^ measurement derived by Becker & Bolton (2013: 
BB13) based on observations of the Lya forest by consid¬ 
ering a wide variety of systematics on this measurement. 
Observations of the Lya forest allow for constraints on 
the emissivity through the impact of the photoionization 
rate r( 0 ) on the prevalence and statistics of Lya-forest 
clouds (and also from constraints on the mean-free path 
Am/p using similar observations). 

We summarize all four of these key constraints in Ta¬ 
ble [T] and Figure [U along with other constraints that 
have been derived in the literature. We refer interested 
readers to R13 and R15 for a comprehensive summary of 
these constraints. 


TABLE 3 

Parameterization for the Cosmic Ionizing Emissivity 
Satisfying the Key Observational Constraints 
Considered Here 



log,Q N\on{z = 8 ) 

(dlog,g Nion/dz)^^s 

Chii 

[s”* Mpc”®] 


2 

50.62«;°5 

_0 1^+008 

3 (fiducial)* 

50.67tO;°| 

a 1C+0.08 

u.io_o 11 

3 (flat prior**)* 

50.651°;“® 

—0 -10+0.09 

5 

50.751°;°® 

—0 14+0-08 

l-|-43z“*-’*’*t 

50.631°;°® 

—0 -l«+0-07 
U.10_o.io 


Excluding Lya Prevalence-Type Constraints'^ 

(i.e., Excluding Constraints 3-4 from Table [1} 

3 50.69t°“« -O.Ogtn® 

More Highly Ionized Universe at 2 : > 7 than S14^ 

3 50.72tg;g^ -0.051^4 

Less Highly Ionized Universe at z 7 than S14^ 

3 50.6lt°;5I -0.20to[\2 

WMAP r = 0.089 zb 0.014, Reionization at z = 5.9-6.5 

3 _ 50.7ltg;g^ _ Q-OOtg-.gg 

* For comparison, we note that only using the Planck 
r = 0.066 zb 0.013 optical depth constraint (PC15), requir¬ 
ing that reionization end at z = 5.9-6.5, and assuming 
that logj^Q iVion( 2 : = 4.75) = 50.99 zb 0.45 (BB13), we find 
logic Mon(z = 8) = 50.66l[(;5g and (d log,Q iVion/cl2)z=8 = 
-O-ieto.ls foi- Chii = 3. 

** In our fiducial determinations, the regions of parameter 
space are weighted according to |Vr| (where the derivatives 
are with respect to log,Q and (d log,g 

Hence, the prior is fiat in units of optical depth r. 

1 Redshift Dependence found in the hydrodynamical simula¬ 
tions of Pawlik et al. (2009). 

As illustrated in Appendix B and Figure [TT] 

Qhii{z = 7) = 0.70 ± 0.15, Qhii{z = 8) = 0.55 ± 0.20. 
Shown in the lower-left panel of Figure [2] with the green- 
shaded region. 

Qhii{z = 7) = 0.49 ± 0.15, Qhii{z = 8) = 0.15 ± 0.20. 
Shown in the lower-left panel of Figure [2] with the orange- 
shaded region. 

until Qhii = 1 when this equation ceases to be valid and 
the ionizing emissivity impacts not only atomic-hydrogen 
gas in the IGM, but also that in higher-density Lya 
clouds. We take Qhii = 0 at z = 25. The recombi¬ 
nation time tree is as follows (KF12): 


3. MODELING THE EVOLUTION OF THE FILLING 
FACTOR OF IONIZED HYDROGEN Qhii 

Here we describe the simple evolutionary models we 
consider for the evolution of the cosmic ionizing emis¬ 
sivity and which we will compare against observational 
constraints on the ionization state of the universe. 

We model the evolution of the cosmic ionizing emissiv¬ 
ity Nioniz) using just two parameters A and B: 

iVio„(z) = 410^(^-8) (1) 


where A = Nion{z = 8) and B = (dlog^o ./Vion/dz)z=8- 
As in previous analyses (e.g., Madau et al. 1999; Bolton 
& Haehnelt 2007; KF12), we follow the evolution of Qhii 
using the relation 


dQnii ( \ Qhii 


( 2 ) 


= 0.88Gyr (i^)" [j^) " ' 

( 3 ) 

where Chii is the clumping factor of ionized hydrogen 
< n^Hii > / < hhii and Tq is the temperature of the 
ionizing hydrogen gas. We adopt a value of 2 x lO'^K 
for the temperature Tq of the ionizing gas to account for 
the heating of the gas that occurs due to the reionization 
process itself (Hui & Haiman 2003). 

Our calculation of the Thomson optical depths t them¬ 
selves also follow familiar expressions from previous anal¬ 
yses (e.g., KF12): 

lTTZ— QHii{z)aTnHQ + 'nY/‘iX) (4) 

do tl[z) 

where ctt is Thomson cross section and X and Y is the 
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Fig. 3.— Allowed evolution ( 68 % confidence) of the cosmic ion¬ 
izing emissivity {region demarcated by the cyan lines, red-shaded 
region, and regions demarcated by the black and green lines) over 
the redshift range 2 = 6 to 2 = ll, as inferred from the three pri¬ 
mary observational constraints considered here (Fan et al. 2006a; 
S14; PC15) and assuming the clumping factor Chii remains fixed 
at 2, 3, 5, and evolves as 1 + 432:“^'^^ (as found by the hydro- 
dynamical simulations of Pawlik et al. 2009), respectively. The 
light-blue-shaded region indicate the constraints from BB13. The 
normalization of the derived ionizing emissivity only shows a weak 
dependence on the assumed value of the clumping factor Chii-, 
changing by just '~^30% for factor-of-2.5 differences in Chii (be., 
Chii = 2 vs. Chii = 5). Strikingly, the evolution inferred for the 
ionizing emissivity, i.e., (dlogj^Q Nion/d^)z=8i is even less sensitive 
to the adopted value for the clumping factor (see also Table 0 . 

primordial mass fraction of hydrogen and helium. Fol¬ 
lowing KF12, we assume that helium is singly-ionized at 
z > 4 (?7 = 1) and doubly-ionized at z < 4 (rj = 2). 

4. RESULTS 

4.1. Quantifying the Redshift Evolution of the Cosmic 
Ionizing Emissivity 

Here we examine what constraints can be set on the 
evolution of the cosmic ionizing emissivity from the key 
observational constraints we consider. We can do this 
without knowledge of the nature of the ionizing sources, 
due to the simplicity of the basic equations that govern 
cosmic reionization, i.e., Eq. m and Nothing in 

these equations requires knowledge of the nature of the 
ionizing sources. 

To derive constraints on the ionizing emissivity, we 
consider a full two-dimensional grid of plausible val¬ 
ues of (1) the normalization of the ionizing emissivity, 
i.e., log^o-Vion( 2 ) 2=8 (from 49.7 to 51.5 s“^ Mpc“^) and 
(2) the dependence of this emissivity on redshift, i.e., 
(d log;,Q fVion/d-2)z=8 (be., from —0.6 to 0.1). For each 
choice of the normalization and redshift dependence of 
the ionizing emissivity {A and B in Eq. [1} , we investi¬ 
gate whether the assumed emissivity would produce a 


reionization history and Thomson optical depth (com¬ 
puted from a reionization history using Eq. a consistent 
with the key observables we consider (i.e., constraints 1-5 
in Table [Hand Figure |T])- 

We begin by considering the model ionizing emissivities 
allowed assuming a redshift-independent clumping factor 
of 3, but then later explore what the impact would be of 
different clumping factors, as well as considering clump¬ 
ing factors that evolve with redshift. Our choice of 3 for 
the fiducial value of the clumping factor is motivated by 
the results of Bolton & Haehnelt (2007) and Pawlik et 
al. (2009: see also Finlator et al. 2012 and Shull et al. 
2012 ). 

To illustrate the impact that various observational con¬ 
straints have on the evolution of the cosmic ionizing emis¬ 
sivities, we start by considering only a subset of the avail¬ 
able constraints. More specifically, we consider the im¬ 
pact of matching both the Thomson optical depths of 
different microwave background missions (best estimate 
and iltr) and requiring that reionization be complete be¬ 
tween z = 5.9 and z = 6.5 (i.e., when Qhii first reaches 
1 using Eq. la, while taking Chii = 3. §2 describes the 
rationale for the Zreion = 5.9-6.5 constraint. The results 
are shown in the top two panels of Figure [2] 

As the top two panels illustrate, there are clear differ¬ 
ences between the implied evolution of ionizing sources 
for WMAP and Planck. Similar to previous models of 
the ionizing emissivity which approximately match the 
WMAP constraints (e.g., Haardt & Madau 2012; KF12), 
the ionizing emissivity we derive shows no (or a slightly 
declining) evolution from z ^ 11. For the Planck con¬ 
straints on T, a range of different models are possible, 
though the ionizing emissivity in all acceptable models 
increases with cosmic time. 

Interestingly enough, the allowed evolution of the cos¬ 
mic ionizing emissivity with redshift shows excellent con¬ 
tinuity with the emissivity inferred at z = 3-5 from the 
Lya forest by Becker & Bolton (2013) [BB13]. This sug¬ 
gests that the evolution we have derived for the ionizing 
emissivity is plausible @ 

In Figure [21 we only present our constraints on the 
ionizing emissivity at redshifts z = 6-12 where changes 
in the filling factor of ionized hydrogen will affect the 
observations we consider. At redshifts lower than z ~ 
6, the IGM is almost entirely ionized, and any change 
in the ionizing emissivity will have little impact on the 
QHiiizfs we consider or the Thomson optical depths. 

More detailed constraints on the evolution can be 
achieved by examining estimates of the filling factor 
of ionizing hydrogen Qhii derived from studies of the 
prevalence of Lya emission in z ~ 7 and z ^ 8 galax¬ 
ies and by requiring that the evolution of the ionizing 
emissivity be consistent (at Icr) with the z = 4.75 con- 

® While the emissivity we derive [extrapolated to z ~ 4] agrees 
with BB13, it agrees less well with that inferred by KF12 (indi¬ 
cated by the hatched-dotted grey region in the top left panel of 
Figure [21l. This is due to a tension between the ionizing emissivity 
results of KF12 and BB13. In noting this tension, readers should 
realize that all such determinations of this emissivity depend quite 
sensitively on various model parameters like the temperature of the 
IGM and opacity to ionizing photons. For BB13’s determination 
of the ionizing emissivity, it was possible to take advantage of some 
new measurements of the IGM temperature over the redshift range 
z = 2.0-4 .8 (Becker et al. 2011), while also including an account 
for cosmological radiative transfer effects. 
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Fig. 4.— (upper) Observational constraints on the ionizing emissivity Nion(z = 8 ) at 2 : ~ 8 and the evolution in this emissivity per unit 
redshift (dlog;i^Q Nion/dz)z=8 (for Chii = 3). 68 % and 95% confidence intervals based on the Planck Thomson optical-depth constraints 
(red and light-red-shaded regions) and combining the prevalence of Lya emission in 2 ~ 7-8 galaxies as found by S14 and continuity with 
the ionizing emissivity at 2 : = 4.75 as derived by Becker Sz Bolton (2013: black lines). The purple lines bracket the allowed parameters 
assuming reionization is completed between 2 : = 5.9 and 2 : = 6.5 (i.e., when Qhii reaches 1 using Eq. n. The blue circle gives the 
equivalent evolution in the ionizing emissivity from the R15 models, (lower) The relative likelihood of different rates of evolution in the 
emissivity, (dlogj^o Mon/^ 2 ) 2 = 8 ? derived by marginalizing over N[on(z = 8 ), along with the maximum-likelihood value and Icr uncertainties 
(red circle and Icr uncertainties). The dotted red line and open circle show the equivalent likelihoods, for a fiat prior in Nioni^ = 3) and 
(dlogiQ N[on/dz)z=8 (rather than in r: §4.1), respectively. Also presented (magenta line) is the equivalent likelihood using the WMAP 
optical depths and assuming reionization finishes between 2 = 5.9 and 2 = 6.5 (as in the upper left panel of Figure O- The expected 
redshift dependence of the ionizing emissivity for galaxies and quasars is also shown with the shaded grey and green regions, respectively. 


straints from BB13. In the lower two panels of Fig- 
ure m the allowed evolution ( 68 % [left+right] and 95% 
[right] confidence) in Nion{z) is shown based on Qhii's. 
derived in S14 (as specified in §2) and also allowing for 
a higher value of Qhii = 0.49 ± 0.15 at z ^ 7 (e.g., 
as in Caruana et al. 2014) and Qhii = 0.55 ± 0.20 
at z ^ 8 . The confidence regions here are derived by 
calculating the probabilities based on constraints 1 and 
3-5 (Table [J) and then marginalizing over the full two 

dimensional logio.Vion(z = 8 ) vs. (dlog^o Mon/dz) 2=8 
space based on the total computed from constraints 
1 and 3-5 (Table [Ij. In marginalizing over the parame¬ 
ter space, regions where Zreion < 5.9, Zreion > 6.5, and 
logiQiVion(^ = 4.75) > 51.44 are excluded. Individual 
regions of parameter space are weighted according to the 
calculated jVrj (where the derivatives are with respect 
to logigiVion and (dlog;io iVion/dz:)z= 8 )- In doing so, we 
weight regions of parameter space according to the im¬ 
pact they have on physical variables like the Thomson 
optical depth r (which approximately varies according 
to the total output of ionizing photons at z 8-15) and 


give less weight to those regions in parameter space (i.e., 
(dlogiQ iVion/d- 2 ) 2=8 < —0.3) which have less impact on 
the ionizing emissivity at z >> 6 . The 68 % and 95% 
confidence intervals presented in the lower right panel of 
Figure [2] are tabulated in Table [2] 

It is clear that if the universe is only 50% ionized by 
z ~ 8 (in significant contrast to the results of S14), it 
would imply a much higher cosmic ionizing emissivity at 
z > 8 . If we assume that the universe is measured to be 
even more ionized at z ^ 7 and z ~ 8 , i.e., Qhii = 0.9, 
as an even more extreme example (this situation might 
arise if use of the Lya fraction as a probe of Qhii is sub¬ 
ject to large systematic errors: e.g., Mesinger et al. 2015: 
but note also Sobacchi & Mesinger 2015 and Keating et 
al. 2014), we derive log;^Q iVion(-z = 8 ) = 50.74+Q g 3 and 

(dlogjQ iVion/dz)z =8 = —O.O 3 I 0 Q 4 . The inferred cosmic 
ionizing emissivity in this case is 0.07 dex higher at z ^ 8 
and 0.31 dex higher at z 10 than our fiducial determi¬ 
nation. 

We found it challenging to reproduce the < 0.3 filling 
factors of ionized hydrogen Qhii at z ~ 7 found by Ono 
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Fig. 5.— 68% and 95% confidence intervals {red and light red-shaded regions) on the UV luminosity density over the magnitude 

interval —23 < Muv,ab < —17 as derived by B15 in specific magnitude intervals {black crosses with Icr errors) and based on a fit to 
the Schechter parameters {black-hatched region). These constraints on the z > 6 luminosity densities are supplemented by the 2 ~ 9 and 
2 : ~ 10 determinations by Ishigaki et al. (2015) and Oesch et al. (2015), respectively. Also shown in the panel are the ionizing emissivities 
we infer here {Chii = 3), offset by a redshift independent conversion factor log^Q /escCion = 24.95. Even though it seems quite clear that 
UV LF extends faintward of —17 mag (e.g., Alavi et al. 2014; Atek et al. 2015; B15), we include this comparison here to illustrate the 
similar evolution observed when minimal extrapolations are employed, {right) Same as left panel, but to a faint-end limit of —13 mag and 
adopting a conversion factor of log^Q fesciion = 24.50. The UV luminosity density integrated down to —13 mag likely evolves more slowly 
with redshift than to —17 mag, based on the steeper shape of the UV LF at high redshift due to an evolution to a (e.g., Bouwens et al. 
2011; McLure et al. 2013) and possibly M^y (e.g., Bowler et al. 2015; B15). The UV luminosity density integrated to —13 mag is also 
more uncertain than integrated to —17 mag due to the greater extrapolation required. In both the left and right panels, the UV luminosity 
density grows at a similar rate to the inferred ionizing emissivities (§4.2). 


et al. (2012) within the context of our simple model for 
the ionizing emissivity, as it implied optical depths of 
<0.046 (in tension with the Planck results at 1.7a) and 
also extrapolated to z = 4.75 implied ionizing emissivi¬ 
ties (i.e., logio .Vionjs”^ Mpc“^]^ 51.5) in excess of that 
measured by BB13|j 

In Figure [21 evolution in the ionizing emissivity is pre¬ 
sented, allowing for different values of the clumping fac¬ 
tor C'hii and also supposing that the clumping factor 
evolves in redshift as 1 -I- as found by Pawlik 

et al. (2009) in sophisticated hydrodynamical simula¬ 
tions. The primary impact of changes in the clumping 
factor is on the overall normalization of ionizing emissiv¬ 
ity Vion, not its evolution with time. This can be clearly 
seen in the maximum likelihood log,QVion(-z = 8) and 
(dlog,Q Vion/d’Z)z=8 values presented in Table [3] for dif¬ 
ferent values for the clumping parameter. Interestingly, 

® The range of different constraints on Qhii based on the preva¬ 
lence of Lya emission in galaxies (often using substantially identi¬ 
cal observations) illustrate the challenges in deriving these Qhii 
factors, as well as their considerable dependence on the simula¬ 
tions used to interpret the available observations (and indeed it is 
clearly non-trivial to adequately capture the many relevant phys¬ 
ical phenomena that enter into these calculations, e.g., growth of 
structure, star formation, radiative transfer of Lya photons, and 
patchy reionization in the same simulation). 


the difference between the derived evolution of the ion¬ 
izing emissivity with Chii = 5 and Chii = 2 is just 
^30%. This is small (e.g., not a factor of 2.5 as might be 
suggested by the ratio of clumping factors) because reion¬ 
ization appears to be a photon-starved process (Bolton 
& Haehnelt 2007). The overall insensitivity of estimates 
of the ionizing emissivity to the clumping factor Chii is 
noteworthy. 

Our 68% and 95% likelihood constraints on the model 
parameters log^^Nioniz = 8) and (dlog^o fVion/d^)z=8 
are also presented in Figure [4] based on optical depths 
measured from Planck, the Qhii^s estimated from the 
prevalence of Lya in z 7-8 galaxies, and requiring that 
the model emissivity extrapolated to z 4.75 matches 
that derived by BB13. A flat prior is assumed in deriving 
these constraints. The purple lines bracket the region 
allowed for reionization to be completed between z = 5.9 
and z = 6.5. 

Marginalizing over all values of Nion{z = 8), we 
Hnd a best estimate (dlog^o fVion/d-z)z=8 of —O.lblg'ii 
{red circle from the lower panel of Figure H]) based on 
the observed constraints. If we adopt a flat prior in 
(dlog,o Nion/dz)z=s and Nic,n{z = 8) (instead of our hdu- 
cial procedure of taking the prior to be flat in r), we find 
a best estimate (dlog,Q fVion/d^:)2=8 of —O.I9toqi {open 
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Fig. 6 .— (left) 68% and 95% confidence intervals {red and light red-shaded regions) on the evolution of the cosmic ionizingemissivity 
(assuming Chii = 3). The derived evolution presented is the same as in the lower right panel of Figure [2] (see also Tabled for these 
constraints in tabular form). Constraints from BB13 are indicated with the light-blue-shaded region. The equivalent trend in Nion derived 
from R15’s model {blue line) is shown here for context and is included among those ionizing emissivity evolution models preferred at 68% 
confidence, {right) 68% confidence intervals on the f/F—luminosity density over the magnitude interval —23 < Muv,ab < —13 as derived 
by B15 in specific redshift intervals {black crosses with la errors). See Tableland Figure|5]for the calculated luminosity densities to other 
faint-end limits. These constraints on the z > 6 luminosity densities are supplemented by the z ~ 9 and z ~ 10 determinations by Ishigaki 
et al. (2015) and Oesch et al. (2015), respectively, extrapolated to the same limiting luminosity (—13 mag). The black-hashed region shows 
the implied evolution of the galaxy f7y-luminosity density, based on our constraints on the ionizing emissivity and assuming that galaxies 
are the source of this emissivity with some redshift-independent conversion factor logj^Q /esc^ion = 24.50. The steep gray line shows the 
redshift dependence one would expect {with approximately the correct normalization [in units of Nian] relative to other quantities presented 
here) for the ionizing emissivity for quasars using the LEDE-model fit from McGreer et al. (2013: see §4.3). The UV luminosity density 
and ionizing emissivity we infer appear to evolve very similarly (§4.2). 


red circle from the lower panel of Figure SJ. 

It is worth emphasizing that the evolution we derive for 
the ionizing emissivity (dlogj^g ^ion/dz)z=8 is only mod¬ 
erately affected by our use of constraints that depend on 
the prevalence of Lya in z ^ 7 and z ~ 8 galaxies (con¬ 
straints 3-4 in Table [T]). Excluding these constraints, we 
find a best value for (d log^g iVion/d-2)z=8 of —O.OOtgqg, 
<0.5tT different from our fiducial determination. The left 
panel of Figure [TT] from Appendix B explicitly illustrates 
the evolution in the emissivity we would infer doing the 
analysis in this manner. 

4.2. Galaxies as the Primary Source of the Cosmic 
Ionizing Emissivity? 

We now consider whether galaxies could plausibly be 
the primary source of the inferred cosmic ionizing emis¬ 
sivity. To explore this question, we want to compare both 
the evolution of the ionizing emissivity and the evolution 
of the galaxy UV-luminosity density. 

To do this, we first examine the observed UV luminos¬ 
ity density from the B15 luminosity function (LF) study 
to observed magnitude limit —17 mag from current ob¬ 
servations (see Figures 15, 18 and 19 from B15). The left 
panel of Figure [S] shows the observed results for the UV 
luminosity density puv from B15 and from other sources 
(e.g., Oesch et al. 2015; Ishigaki et al. 2015). 


It is clear given the steep slope at —17 mag and the lack 
of any indication of a cut-off from Figures 15 and 19 in 
B15 for example (see also Schenker et al. 2013; McLure et 
al. 2013; Alavi et al. 2014; Barone-Nugent et al. 2015), 
as well as early results from magnified sources found in 
leasing clusters (e.g., Atek et al 2015) that fainter galax¬ 
ies must contribute substantially to the total UV radia¬ 
tion from galaxies (e.g., Yan & Windhorst 2004; Beck¬ 
with et al. 2006; Bouwens et al. 2007; Reddy & Steidel 
2009). While we can make plausible extrapolations based 
on tbe derived Schechter parameters, the question arises 
as to the appropriate limit. 

As others have done (e.g., R13, R15) we adopt a fidu¬ 
cial luminosity of —13 mag down to which galaxies are 
typically assumed to be able to form efficiently. Faint- 
ward of —13 mag, galaxy formation may be suppressed 
due to the inefficient gas cooling onto low-mass halos 
(e.g., Rees & Ostriker 1977) or due to the difficulties in 
low-mass galaxies retaining their gas (e.g., Mac Low & 
Ferrara 1999; Dijkstra et al. 2004). The choice of the 
limiting luminosity to adopt for the LF, e.g., Mum, is 
also motivated from sophisticated hydrodynamical sim¬ 
ulations (e.g., O’Shea et al. 2015) or from fits to the 
faintest points in the LF (e.g., Munoz & Loeb 2011; 
Barone-Nugent et al. 2015). Some guidance can also 
be obtained by attempting to match the redshift de¬ 
pendence inferred for ionizing emissivity Nio^{z) with 
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Fig. 7. — 68% and 95% confidence intervals on parameters describing the cosmic ionizing emissivity (assuming Chii = 3: red contours) 
and a comparison to the emissivity expected (gray boxes in the upper middle) for galaxies to three different faint-end cut-offs to the LF 
(—10, —13, and —16 mag) where the Lyman-continuum photon production efficiency ^jon and fesc plausibly have values over the wide 
range: to s“^/(ergs s“^Hz“^) and 0.05 to 0.40 based on the observations (e.g., Siana et al. 2010, 2015; Vanzella et al. 2012; 

Mostardi et al. 2013; Cooke et al. 2014; Bouwens et al. 2014). The arrows indicate how potential evolution in the C/l/-continuum slopes 
^’s or fesc of galaxies (as d^jdz —0.04 [predicted by Wilkins et al. 2013 and consistent with the evolution observed by Bouwens et al. 
2014] or fesc oc (1 -t- required by HM12 to match the WMAP r measurements) would affect the evolution of iVion- The gray box to 

the lower left shows the expectations for quasars using the Willott et al. (2010) results to set the normalization of the ionizing emissivity 
N[on(^) (after correcting their results upwards by 0.2 dex to account for a possibly steeper faint-end slope a ~ —2 versus the a = —1.5 
Willott et al. 2010 adopt) and the redshift dependence of the emissivity from the McGreer et al. (2013) LEDE and modified LEDE models. 
The 68% likelihood contours preferred based on the WMAP opticle depths and assuming reionization finishes a.t z = 5.9-6.5 are shown 
in magenta. The small size of the WMAP contours is an artifact of the sensitivity of the computed optical depths r to small changes in 
(dlog^O-^ion/c^- 2 :) 2=8 when (dlogiQ N[on/dz)z=8 approaches 0 (since such implies a constant puv to arbitrarily high redshift). The gray 
and blue circles give the equivalent parameters from the HM12 and R15 models, respectively. It is clear that galaxies should be successful 
at producing the inferred ionizing emissivity for a variety of plausible values for /esc, ^ion, ^-i^d Mn^. It is also clear from this figure that 
quasars/AGN seem quite unlikely to be the source of this emissivity (see also HM12). 


Puv (e.g., Bouwens et al. 2012a; R13), i.e., compare 
the inferred (dlog 2 Q_yion/d 2 :)z= 8 ’s in Table [3] with the 
dpuv/dz's in Table SI but this presupposes that galaxies 
are the only source of the ionizing UV radiation and so 
it is considered indicative. 

With our adopted limiting magnitude of —13 mag, we 
estimate the UV-luminosity density puv implied by the 
B15 LF results, by marginalizing over the full likelihood 
distribution in M*, (jf, and a and computing both the 
mean and la error from the resultant likelihood distri¬ 
bution for pi/v (similar to Figure 3 from R13). 

The estimated UV-luminosity densities puv to —13 
mag are shown in the right panel of Figure [5l The shal¬ 
lower redshift dependence of the result is due to the inclu¬ 
sion of much lower-luminosity galaxies in the calculated 
luminosity densities and the strong evidence that the UV 
luminosity density evolves faster at the bright end than 


the faint end (e.g., Yoshida et al. 2006; Bouwens et al. 
2006, 2008; B15). The most important factor in this dif¬ 
ferential evolution is the flattening of the faint-end slope 
a to the UV LF with cosmic time (e.g., Bouwens et al. 
2011; Schenker et al. 2013; McLure et al. 2013). For con¬ 
text, the luminosity densities puv implied by the B15 
Schechter parameters to many different limiting lumi¬ 
nosities and as a function of redshift are presented in 
Table ID The red regions on the two panels of Figure [5] 
will be discussed below following the discussion of Fig¬ 
ure El 

We are now positioned to explore both the evolu¬ 
tion of the ionizing emissivity and the evolution of the 
galaxy UV-luminosity density, as shown in the left and 
right panels of Figure El respectively, with the vertical 
axes chosen so that the inferred ionizing emissivity and 
luminosity density maximally overlap, i.e., iVion(2) = 
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As is apparent from Figure [ 6 l the best-fit evolu¬ 
tion in the ionizing emissivity (d logj^g-^ion/rf-z)z =8 of 
—0.15 ^q’5i (or —0.19;l^g °® adopting a flat prior in 

(dlogio-/Vion/(i 2 )z =8 and logioiVion(z = 8 )) is in excel¬ 
lent agreement with the best-fit evolution in the UV 
luminosity density to —13 mag, i.e., dlogiQ puv/dz = 
—0.11 ±0.04. While the uncertainties are still large, this 
is suggestive that early star-forming galaxies provide the 
ionizing photons needed to reionize the universe. 

We emphasize that the present conclusions do not sig¬ 
nificantly depend on our use of the observational con¬ 
straints based on the prevalence of Lya emission from 
z ~ 7 and z ~ 8 galaxies. As we discuss in §4.1, if we 
exclude these constraints, the (dlog^g 77ion/d2:)z=8 value 
we derive (—O.OOt^ia) almost identical to the best- 
fit evolution we find in the UV luminosity density. See 
Appendix B and Fig. [TT] where this point is illustrated 
explicitly. 

It is interesting also to consider the offset where the 
ionizing emissivity and UV luminosity density maximally 
overlap. The best overlap occurs adopting an offset of 
1024.50 s-i/(ergs s“^Hz“^) (= /esc 6 on) for a faint-end 
limiting luminosity of —13 mag (similar to what has been 
proposed in many other studies: e.g., KF12, R13). This 
conversion factor has an uncertainty of at least 0.10 dex, 
given the uncertainties on both the normalization of the 
ionizing emissivity at z 8 (typically ^0.08 dex: Ta¬ 
ble [3]) and the ±0.05 dex uncertainties on the UV lumi¬ 
nosity density integrated to < —13 mag (see Table a. 
This is one of the first analyses to estimate an uncer¬ 
tainty on this conversion factor using observations that 
directly concern the z > 6 universe (see also Mitra et al. 
2013). The evolution of the ionizing emissivity is also 
shown in relation to the puv for the two magnitude lim¬ 
its in Figure [Sj using an appropriate offset. It is quite 
clear that the two quantities seem to evolve in a similar 
manner even to the bright limit. 

Is this multiplicative offset plausible if star-forming 
galaxies are to be the source of the cosmic ionizing emis¬ 
sivity? Given that it can be expressed as the product 
of two factors ^ion, the production efficiency of Lyman- 
continuum photons per unit UV luminosity, and the es¬ 
cape fraction /esc with plausible values in the range 10 ^^'^ 
to 10^®-^ s“^/(ergs s“^Hz“^) and 0.05 to 0.40 based on 
observations of z ^ 2-4 galaxies (Vanzella et al. 2012; 
Mostardi et al. 2013; R13; Duncan & Conselice 2015), the 
observed offset between the inferred ionizing emissivity 
and UV luminosity density is certainly within the range 
allowed by these value, i.e., 10 ^^ ® to 10 ^^-^ s“^/(ergs 
s“^Hz“^). This is the case for all three limiting lumi¬ 
nosities Miim presented in Figure[7](— 10 mag, —13 mag, 
and —16 mag). 

In distinct contrast to the situation faced in interpret¬ 
ing the WMAP r’s, no evolution in the escape fraction 
is required to match the ionizing emissivity model we in¬ 
fer, as had been earlier considered by HM12 or KFI2. 
Neither is evolution in the C/U-continuum slopes /3 re¬ 
quired, despite the apparent trend towards bluer /3’s at 
high redshift (e.g., dfdjdz ^ —0.10 ± 0.05: Bouwens 
et al. 2014: see also Bouwens et al. 2012; Kurczyn- 
ski et al. 2014; Finkelstein et al. 2012a; Wilkins et al. 


TABLE 4 

XJV Luminosity Density derived from the B15 LF 

PARAMETERS INTEGRATED TO DIFFERENT FAINT-END 
CUT-OFFS ALjim TO THE UV LF. 


^lim 

logio 

Puv [ergf 

i s-iRz' 

^ Mpc 


-10 

(26.47 

± 0.08) - 

(0.07 

± 

0.05)(z - 

-6) 

— 13 (fiducial) 

(26.37 

± 0.05) - 

(0.11 

± 

0.04)(z - 

-6) 

-16 

(26.20 

± 0.03) - 

(0.16 

± 

0.03)(z - 

-6) 

_17a 

(26.11 

± 0.03) - 

(0.19 

± 

0 .02)(z - 

-6) 

-17.7=' 

(26.02 

± 0.03) - 

(0.21 

± 

0 .02)(z - 

-6) 


The evolution of the UV luminosity density to —17 
and —17.7 mag (as used by Bouwens et al. 2011; Ellis 
et al. 2013; Oesch et al. 2013, 2014, 2015; McLeod et al. 

2015) is included for contrast with the evolution seen 
to much fainter limits and also for completeness. The 
steep slope of the UV LF at the —17 mag limit of deep 
searches for faint galaxies at z ~ 7-8 suggests that the 
cut-off is fainter than —17 mag (e.g., B15; Schenker 
et al. 2013; McLure et al. 2013). Recent results from 
z ~ 7 and z ~ 8 LFs derived from magnified sources 
behind lensing clusters also indicate a steep slope, and 
extends the limit at z 7 to —15.5 mag (e.g., Atek et 
al. 2015). 

2011 ) or as expected from simple theoretical models (e.g., 
d/3/dz —0.04 [Wilkins et al. 2013] or d/3/dz ^ —0.1 
[Finlator et al. 2011]). A mild evolution in ^ is not in¬ 
consistent with our constraints on the ionizing emissivity. 

Given the likely connection between the ionizing emis- 
sivities we infer and galaxy UV luminosity density at 
z > 6, what are physically plausible values for ^lon and 
/esc that we would infer? The relevant value of the C/U- 
continuum slopes /3 of galaxies to estimate their ionizing 
emissivity contribution is ^ —2.3, as most of the luminos¬ 
ity density at z > 6 occurs in lower luminosity galaxies 
and the median /3 for faint z 7-8 galaxies is —2.3 
(e.g., Bouwens et al. 2014). Using a similar approach 
to R13, one can show that /ion in such galaxies has an 
approximate value of 10 ^^-"^® s“^/(ergs s“^Hz“^) for /ion 
(Appendix A), similar to the /ion advocated by Topping 
& Shull (2015). For this value of /ion, the escape frac¬ 
tion fesc is 0 . 1 llo;o 2 for the observed conversion factor 
1024.50 g-i/^ergs s“^Hz“^). Alternatively, if we take /ion 
to be equal 10^®'^ s“^/(ergs s“^Hz“^) as R13 adopted to 
match with the /3 measurements of Dunlop et al. (2013), 
the relevant value of fesc is 0 . 20 )to)o 4 (essentially identical 
to R13’s adopted value of 0.2). 

Glearly, many degenerate combinations of fesc, /ion, 
and Miim can be successful in producing the same ioniz¬ 
ing emissivity. For variable Mum, the present constraint 
on /esc/ion Can be generalized to the following constraint 
on these three parameters: 

fescfionfcorriMiim) = 10^^®°S“V(CTgS S“^Hz“^) (5) 

where the added term fcorr{Miim) = 
2Q0.02-i-0.078(MHm-i-i3)-0.0088(MHm-i-i3)^ corrects 

Puv{z = 8 ) derived to a faint-end limit of Mum = —13 
mag to account for different Mum's 0 For clumping 
factors Chii not equal to our fiducial value of 3, the 

In deriving this correction factor, we made use of the following 
relationship between puy {z = 8) and the faint-end cut-off to 

the LF: puv{z = 8) = (26.17 ±0.09) ± (0.08 ± 0.02)(MHm ± 

13) — (0.009±0.008)(M/i^ ± 13)^. This relationship can be derived 
by fitting to the results in Table [4| and is accurate to 5%. 
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left-hand side of the above equation should be multipled 
by based on the results presented in 

Table[3](the {Chii scaling is approximately valid 
for Chii < 10). The factor s“^/(ergs s“^Hz“^) 

has an uncertainty of ~0.1 dex. 

Eq. m can be manipulated to allow for an estimate 
of fesc given assumed values for the other parameters 
Miim, ^ion, and Chi I, assuming that galaxies reionize 
the universe: 

fesc « C^/-J,(M;,„)(Cff„/3)°-3l024-5°s-V(ergs s-^Hz-i) 

( 6 ) 

Estimates of fesc for various fiducial choices of ^ion, 
Miim, and clumping factors Chii are provided in Ta¬ 
ble [5l This is one of the first analyses to estimate the 
uncertainty on the derived value of fesc based on obser¬ 
vational constraints on the ionization state of the z > 6 
universe (see also Mitra et al. 2013). 

The issue of degeneracies among the parameters /esc, 
^ion, and Miim is discussed extensively in KF12 (see also 
Bouwens et al. 2012a; Finkelstein et al. 2012b; R13) and 
may not be easy to resolve based on observations in the 
immediate future. 

4.3. Quasars as the Primary Source of the Cosmic 
Ionizing Emissivity? 

We briefly consider whether quasars could be the pri¬ 
mary source of the cosmic ionizing emissivity. Despite 
their relative scarcity in the z > 4 universe, quasars 
potentially can contribute quite substantially to the in¬ 
ferred emissivity due to the hardness of their spectrum 
and a much higher escape fraction {fesc = 1?: e.g., Loeb 
& Barkana 2001). 

While a fraction of the emissivity from quasars at z > 4 
is expected to originate from the most luminous sources, 
a potentially large fraction of their contribution could 
originate at much lower luminosities, and therefore it is 
important to have reasonable constraints on both the 
volume densities of faint quasars (including AGN) and 
their faint-end slopes to effectively estimate the ionizing 
emissivity they produce. 

Perhaps the deepest, wide-area probes of the z > 4 
quasar LFs are provided by Willott et al. (2009, 2010) 
and McGreer et al. (2013). The deepest part of Willott 
et al. (2009, 2010) searches for faint (> —22 mag) z ~ 6 
QSOs over 4.5 deg^, while McGreer et al. (2013) lever¬ 
aged the deep observations over a 235 deg^ region in 
SDSS Stripe 82 to probe the prevalence of moderately 
faint (> —24 mag) z ^ 5 qusars. Willott et al. (2009) 
identified 1 very faint z = 6.01 quasar over the CFHTLS 
Deep/SXDS, while McGreer et al. (2013) identified some 
^70 faint quasars in their search, 29 with absolute mag¬ 
nitudes faintward of —25 mag. 

Given the relatively small numbers of faint quasars 
identified by these programs (and other recent searches 
e.g., Weigel et al. 2015, who find no convincing z > 5 
AGN over the Chandra Deep Field South), the total 
emissivity of ionizing photons from the quasar popula¬ 
tion at z ^ 4-6 is still quite uncertain. Nevertheless, a 
sufficient number of faint quasars have been found that 
estimates of the total emissivity of the population can 
still be made (but see the discussion in Giallongo et al. 
2015). At this time, the most reliable estimates can 
be made using the faint-end slopes a and evolutionary 


trends derived by McGreer et al. (2013). McGreer et al. 
(2013) find a log^g 4>{z) = log^g (jio — 0.60(z — 2.2) and 
M^,{z) = Mtffi — 0.68(z — 2.2) trend with their LEDE 
model and a log^g (j){z) — log^g (j)o — 0.70(z — 2.2) and 
M^{z) = M*^g — 0.55(z — 2.2) trend with their mod¬ 
ified LEDE model. These models imply a comoving 
emissivity log^g e(z) that evolves as —0.33(z — 2.2) and 
-0.48(z-2.2). 

If we take the ionizing emissivity estimates that Willott 
et al. ( 2010 ) derive from their deep z ^ 6 search, extrap¬ 
olate to z ~ 8 , and adjust their result upwards by 0.2 
dex to account for the steeper faint-end slope of a = —2 
found at z ^ 5 by McGreer et al. (2013: versus an as¬ 
sumed faint-end slope of o = —1.5 by Willott et al. 2010), 
the logjQ fVion( 2 : = 8 ) we would derive is 48.5, with a 
(d logjg iVion/dz) 2=8 of —0.41±0.08 taking the arithmetic 
mean of the two LEDE models considered by McGreer 
et al. (2013). Fan et al. (2001) derive —0.47±0.15 based 
on the evolution of the space density of bright quasars 
from z = 6 to z = 3. If we compare these parameters 
with the values we derive for the ionizing emissivity, it is 
clear that faint quasars do not appear to come close to 
providing enough photons to drive the reionization of the 
universe. In addition, the ionizing emissivity produced 
by this population shows a redshift dependence that is 
slightly steeper than what we require to reproduce the 
observational constraints (see e.g. the green-shaded re¬ 
gion and the red line at the bottom of Figure 0]). 

Is it possible that even fainter quasars/AGN could 
boost these luminosity densities? Given that the faint- 
end slope a McGreer et al. (2013) find at z ~ 5 is 
—2.03/g :[ 4 , the luminosity density will be only logarith¬ 
mically divergent, which for an integration to ~ —22 mag 
and ^ —16 mag (if such fainter quasars/AGN exist in 
large numbers) would only increase the overall luminos¬ 
ity density by a factor of 1.5 (0.2 dex) and a factor of 4 
(0.6 dex). This suggests that a full consideration of the 
contribution from faint quasars can potentially boost the 
total ionizing emissivity produced by quasars. However, 
even with such steps, quasars appear quite unlikely to 
contribute meaningfully to the reionization of the uni¬ 
verse at z > 6 (but see Giallongo et al. 2015). 

4.4. Implications for the UV Luminosity Density at 
z - 10? 

The current discussion and others (e.g., R15; Mitra et 
al. 2015) suggest that star-forming galaxies produce the 
bulk of the ionizing emissivity. Given this, we can use 
current constraints on the evolution of the ionizing emis¬ 
sivity to provide an estimate of the UV luminosity den¬ 
sity at z ~ 10. This is useful since there has been much 
discussion about whether the UV luminosity and SFR 
density are trending differently with time from z ^ 10 to 
z ~ 8 than at later times (e.g., Oesch et al. 2012, 2013, 
2014, 2015; Ellis et al. 2013; Goe et al. 2015; Ishigaki et 
al. 2015; McLeod et al. 2015). 

In the simplest case of no evolution in fesc or ^ion, evo¬ 
lution in the UV luminosity density would mirror that 
seen in ionizing emissivity. In this case, one would ex¬ 
pect a factor of increase in the UV 

luminosity density from z ^ 10 to z ^ 6 integrated to the 
faint-end cut-off to the LF Mum, which we will take to 
be —13 mag. However, since the UV luminosity density 


TABLE 5 

Required Values of /esc for different Mh^, ^ion, and clumping factors Chii assuming that 

GALAXIES drive THE REIONIZATION OF THE UNIVERSE^ 
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^ion — 10 

■is.ib s-J-/(ergs 

s-lRz-l) 

^ion — 10 

■ib.-z s-J-/(ergs s’ 


Chii 

Mlim = —17 

Mlim = -13 

3 

II 

1 

0 

Mlim = -17 

Mlim = —13 

Mlim = —10 

2 

u.^o_o,o5 

'^•^'^- 0.02 

0 

'J-O»o_o.oi 


n 1 q+0.05 
'1--Lo_0.04 

^•-*-^- 0.02 

3 

u.^y _0 06 

0 11+0.03 

^•-LJ-_o.02 

O'O+o.M 

0.52l«;)4 

^•^^-0.04 

0 13+0.03 

5 

0 34+0 09 

0 13 + 0.03 

0.08l°°2 

0'-OJ-_o.i3 

n 90+0.06 

'^•^'J_0.05 

n l^+0.04 


-0.06 

0 10+003 
^•^^- 0.02 

O-O+S.o? 

^•^y o.io 

0 19+005 
^•^^-0.04 

^•-^^-0.03 


^ Redshift Dependence found in the hydrodynamical simulations of Pawlik et al. (2009). 

^ These fesc factors can be derived from Eq. Ill in §4.2 of this paper. Importantly, we can also quote 
uncertainties on the estimated /esc’s, which follow from our Icr error estimate (~0.1 dex) on the conversion 
factor 10^'^'^®s“^/(ergs s“^Hz“^) from UV luminosity density puy to the equivalent ionizing emissivity 
Mon- Constraints on fesc are also attempted by KF12 at 2 : 4 based on the derived Mon there (see also 

Finkelstein et al. 2012b for an estimated fesc at 2 ~ 6 based on Mon from Bolton & Haehnelt 2007). 

^ Adopted by R13. 

^ In fact, this is the same fesc that R13 and R15 suggest using for the fiducial parameter choices they 
adopt: ^ion = 10^^’^ s“^/(ergs s“^Hz“^), Mn^ = —13, and Chii = 3. However, no attempt was made 
by R15 to provide a constraint on the uncertainties in this estimate of fesc as derived here. 
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Fig. 8.— Constraints on the UV luminosity density at 2 10 to 

the typical observed limit of —17 mag {red circle), derived from the 
2 = 6-10 evolution of the ionizing emissivity (§4.4). The luminosity 
density results from B15 at 2 = 6 are used as a baseline. The 
derived UV luminosity density at 2 = 10 is also corrected downward 
by —0.32 dex to account for the different evolution in puy expected 
to —17 mag vs. to —13 mag (see Table[^. For context, estimates of 
the UV luminosity density to the same faint-end limit based on the 
Reddy ^ Steidel (2009) 2 = 3 results {green cross), the B15 2 = 4- 
10 results {blue circles), the McLure et al. (2013) 2 = 7-8 results 
{cyan circles), the Ishigaki et al. (2015) 2 = 9 {cyan open circles), 
and Oesch et al. (2015) 2 = 10 results {blue open circles) are also 
shown. The shaded-light-blue curve indicates the evolution in UV 
luminosity density one would expect based on extrapolations of 
the B15 LF parameters to 2 > 8. The dotted green Icr upper limit 
at 2 = 11.5 shows constraints on the UV luminosity density at 
2 > 10 from R15 derived using the recent Planck r measurements 
corrected to —17 mag based on the empirical scalings presented in 
Table El 

is measured directly from observations to about —17 AB 
mag in the reionization epoch, it is perhaps preferable to 
evaluate the changes to that well-determined limit and 
to correct the derived results from the ionizing emissivity 
to —17 mag instead of the extrapolated —13 mag limit. 
Considering the evolution integrated to —17 mag, we 
would expect dlog,Q pc/y= —0.19 — (—0.11) = 0.08 


more evolution at the bright end than the faint-end (us¬ 
ing the scalings from TablelD see also Yoshida et al. 2006; 
Bouwens et al. 2007, 2008 where such differential trends 
were first discussed). This would suggest an increase of 
81*14® X the luminosity density from z ~ 10 to z ^ 6 
to —17 mag. Integrating the z ^ 6 B15 LF to —17 mag 
and accounting for this evolution, we estimate a UV lu- 
minosity density of 10 ^-0.44 ergs s“ Hz“ Mpc“'^ at 

z ^ 10 to —17mag (Figure [8|). This estimate should be 
regarded as an upper limit, as any expected mild evolu¬ 
tion in fesc or ^ion (Siana et al. 2010; Hayes et al. 2011; 
Bouwens et al. 2014: but see also R13) would likely act 
to lower the derived luminosity density at z ^ 10. 

How does this luminosity density compare with a sim¬ 
ple extrapolation of the z = 4-8 LF results to z ^ 10? 
Adopting the dlog4o Pc/y/dz: = —0.19 ± 0.04 scaling im¬ 
plied by the LF results of B15 (Table|4]), the extrapolated 
LF density at z ^ 10 is gj.gg g-i Hz“^Mpc“® 

to —17mag (see light-blue-shaded contour in Figure [5]). 
The luminosity density we infer is consistent with this ex¬ 
trapolation (similar to recent results by Coe et al. 2013, 
McLeod et al. 2015, or the Oesch et al. 2015 results over 
the first Frontier Field). However, it is also consistent 
at 1(7 with the (1 -I- z)“^° ® evolution found by Oesch et 
al. (2014) at z > 8, which suggests a luminosity density 
2 q25.i±o. 3 gj.gg g-i Hz“^Mpc“®. This is particularly the 
case, since z > 6 galaxies may be more efficient at re¬ 
leasing ionizing radiation into the IGM due to evolution 
in fesc or ^ion (both of which have been speculated to 
increase at z > 6: Siana et al. 2010, Hayes et al. 2011, 
HM12, KF12, Duncan & Conselice 2015). This makes 
the present estimate of the UV luminosity density at 
z ~ 10 effectively an upper limit. 

4.5. How the Ionizing Emissivities We Infer Compare 
with Key Observational Constraints 

Finally, it is useful to compare the results of our pre¬ 
ferred models for the evolution of the ionizing emissivity 
with the key observational constraints we considered, as 
a check on the overall self-consistency of the constraints. 

The results for our models preferred at 68% and 95% 
confidence and where reionization finishes at z = 5.9-6.5 
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Fig. 9.— Comparison of the key observational constraints considered here with the results from the simple two-parameters models for 
the cosmic ionizing emissivity preferred at 68% and 95% confidence (§4.5). (left) Shown are the constraints on the Thomson optical depth 
T provided by the Planck 3-year results (PC15: cross-hatched black region). The red and light-red-shaded regions show the range of 
cumulative Thomson optical depths for our models of the ionizing emissivity preferred at 68% and 95% confidence, respectively (Figure [4jl, 
and where reionization is complete between 2 : = 5.9-6.5. (right) Shown are constraints on the filling factor of ionized hydrogen Qhii 
as a function of redshift. The constraints are largely as compiled by R15 (see Tabled of the present manuscript) and are based on the 
Gunn-Peterson optical depths and dark-gap statistics measured in 2 : ~ 6 quasars (Fan et al. 2006a; McGreer et al. 2015: solid black circles 
and squares), damping wings measured in 2 ~ 6.2-6.4 quasars (Schroeder et al. 2013: open black square) and a 2 = 6.3 GRB (Totani et al. 
2006; McQuinn et al. 2008: open green circle), Lya-emitter LFs and clustering statistics at 2 6.6 (Ouchi et al. 2010: open blue square 

and circle) and at 2 7 (Ota et al. 2008: open blue circle), and the prevalence of Lyo emission in 2 ~ 7-8 galaxies (S14: large red circles). 

Also included here (small red solid circles) are other estimated constraints on Qhii from the prevalence of Lyo emission from galaxies at 
(Qhii < 0.5 [R13]; Qhii < 0.4 [Ono et al. 2012]; Qhii < 0.49 [P14]; Qhii ^ 0.5 [Garuana et al. 2014]) and at 2 8 (Qhii < 0.7 

[Tilvi et al. 2014]). The red and light-red-shaded region indicates the range of Qhii allowed for our models of the ionizing emissivity 
preferred at 68% and 95% confidence, respectively, and where reionization is complete between 2 = 5.9 and 2 = 6.5. The magenta-hatched 
region indicates the range of Qhii allowed at 68% confidence for the WMAP 9-year r measurement (0.089 ib 0.014). 


(i.e., when Qhii reaches 1 using Eq. [2]) are presented in 
Figure[9] Both the preferred optical depths and reioniza¬ 
tion histories are shown in this figure. Our ionizing emis¬ 
sivity models are in excellent agreement with all available 
constraints, including those not considered in deriving 
the ionizing emissivity (Table [T]). A similar version of 
this figure is presented in Appendix C, but also showing 
the constraints based on the WMAP 9-year r’s. 

As in other simple models of cosmic reionization (R13, 
R15), i.e., where inhomogeneities in the IGM are not 
treated, we are not able to reproduce observations which 
suggest incomplete reionization to z = 5-6, i.e., Qhii ~ 
1 — 10“^. Correctly reproducing the end of reionization 
would require careful simulations over cosmologically- 
significant volumes, with voids and overdensities, and 
require the treatment of radiative transfer effects. 

Interestingly enough, the optical depths allowed at 68% 
confidence only include the lower 50% of the values al¬ 
lowed by Planck and do not significantly extend above 
r’s of 0.066. This appears to be the direct result of the 
constraint we apply on the filling factor of ionized hydro¬ 
gen at z = 8 from S14. 

If we do not consider the constraint from SI4 and in¬ 
stead suppose that Qhii = 0.55±0.20 at z ^ 8 (resulting 
in an evolution of the ionizing emissivity illustrated by 


the shaded-green region in the lower-left panel of Fig¬ 
ure [I, the Thomson optical depths we derive span the 
range 0.060 to 0.082. This demonstrates the value of con¬ 
tinued work on quantifying the prevalence of Lya emis¬ 
sion in z ~ 8 galaxies (e.g., S14) and also constraining 
the ionizing state of the universe at z ~ 8 using other 
methods (e.g., as Bolton et al. 2011 do at z ~ 7.1 using 
the damping wing of Lya observed against the spectrum 
of a bright quasar). 

5. DISCUSSION AND SUMMARY 

Here we have combined the new measurements of the 
Thompson optical depth r from Planck (PC15) with a 
new approach that focuses on inferring the cosmic ioniz¬ 
ing emissivity Mon and its evolution over the reionization 
epoch from a variety of observables (building of course 
on significant earlier work by e.g. Madau et al. 1999; 
Bolton & Haehnelt 2007; KF12; R13; Becker & Bolton 
2013; R15: see also Mitra et al. 2011, 2012). This ap¬ 
proach has allowed us to first gain insight into the allowed 
evolution of the sources driving reionization, without im¬ 
mediately making assumptions about their nature. We 
then assess the implications of the derived evolution for 
the ionizing emissivity and compare with the evolution 
seen in the luminosity density puv for galaxies. 
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We have derived empirical constraints on the evolution 
of the ionizing emissivity based on the recently-measured 
Thomson optical depth r (PC15) and the QuiiizYs in¬ 
ferred from quasar absorption spectra and the prevalence 
of Lya emission in z = 7-8 galaxies (Fan et al. 2006a; 
McGreer et al. 2015; S14). We tabulate in Table [T] the 
extensive results on the filling factor Qhii of ionized hy¬ 
drogen from the literature that allow for us to constrain 
the evolution of the ionizing emissivity iVion (see also 
Figure [J). 

We demonstrate that the evolution in the cosmic ion¬ 
izing emissivity at z > 6, i.e., (d logj^Q Mon/d-z)z=8 = 
-0.15tn? ((dlogioMo„/dz),=8 = -0.19;°:™ for a 
flat prior in (d log^o Mon/d2:)2=8 and logioMon(2 = 
8)), is matched by similar evolution in the UV lumi¬ 
nosity density (after extrapolation to —13 mag), i.e., 
AlogiQ puv/dz = —0.11 ± 0.04. By contrast, use of the 
9-year WMAP optical depths r = 0.089 ± 0.014 (Ben¬ 
nett et al. 2013) to derive the evolution in the emissivity 

yields (dlog^gfVion/d^)2=8 = 0.00to;o6 (TableE]). 

This is the first time this similar evolution has been 
shown quantitatively in this manner and builds on the 
well-known case that galaxies taken together can provide 
the UV ionizing radiation needed for reionization (but 
see also R15). This result further supports the view that 
star formation in early galaxies drives the reionization of 
the universe (§4.2). 

This conclusion is further strengthened by the fact that 
the requisite conversion factor from the UV luminos¬ 
ity density to the ionizing emissivity (Cion/esc = 10^^'^° 
s“^/(ergs s“^Hz“^)) is consistent with physically plausi¬ 
ble values for the escape fraction fesc and ^ion, for a faint- 
end limit Miim to the UV LF of —13 mag and a clumping 
factor Chii = 3. We calculate that this conversion factor 
has an approximate uncertainty of ^0.1 dex (§4.2) based 
on the uncertainties in the normalization of the cosmic 
ionizing emissivity and the UV luminosity density. We 
also present a generalization of this constraint for other 
values of Mum and Chii, be., Eq. ®. We manipulate 
this constraint to provide a general formula for deriving 
fesc for a wide range of different values for f,ion, Mum, 
Chii assuming that galaxies drive the reionization of the 
universe (see Eq. Inland Table [5]). 

This is also one of the first analyses where the un¬ 
certainty on this conversion factor has been estimated 
from constraints on the ionization state of the universe 
at z > 5.9 and is the direct result of the inferences we 
make regarding the evolution of the ionizing emissivity 
(but see also Mitra et al. 2013). [KF12 achieved a sim¬ 
ilar constraint on the conversion factor at z ~ 4 from 
the ionizing emissivity estimates based on observations 
of the Lya forest.] Despite the size of these uncertain¬ 
ties, the consistency with physically plausible values of 
fesc, 6on, and Miim was not assured a priori and pro¬ 
vides confidence that galaxies play a dominant role in 
reionization. 

We also consider quasars/AGNs as potential sources 


of the ionizing radiation. However, as for most previ¬ 
ous assessments, quasars/AGNs, appear unlikely to be 
the dominant source of the ionizing UV radiation, even 
under rather generous assumptions about the contribu¬ 
tions of faint AGNs. There is simply little evidence they 
show the required emissivities nor redshift dependence to 
match that found for the ionizing emissivity (§4.3: but 
see Giallongo et al. 2015). 

Assuming no change in the production efficiency of UV 
continuum photons or ionizing photons, our constraints 
on the evolution of the ionizing emissivity can be used 
to estimate the UV luminosity density at z ^ 10. We 
show in §4.4 that it is 8^4^ x lower than at z ^ 6. As we 
note, there is also the possibility, due to either evolution 
in Cion or fesc, that z > 6 galaxies are more efficient 
(per unit UV luminosity) at releasing Lyman-continuum 
radiation into the IGM than z 6 galaxies. If that is the 
case, then this estimate on the UV luminosity density of 
z ov/ 10 galaxies is an upper limit. The uncertainty on 
this estimate, however, is large enough that it cannot 
help to resolve the question about the potential change 
in the slope of the luminosity density at around z ^ 8 as 
first identified by Oesch et al. (2012). 

The recent remarkable observations from Planck 
(PG15) provide a fresh opportunity to re-evaluate the 
role of galaxies in cosmic reionization (see also R15; Mi¬ 
tra et al. 2015). We have taken advantage of these new 
Planck results on the Thompson optical depth, as well 
as a decade of observations of galaxies and quasars that 
provide constraints on the filling factor of ionized hydro¬ 
gen as a function of redshift (see §4.5), to take a further 
step. We show not only that reionization is consistent 
with being driven by the UV radiation from galaxies, as 
many others have demonstrated, but also that the evo¬ 
lution of the ionizing emissivity (from z ~ 4 to z ~ 10) 
matches a similar trend in the UV luminosity density 
(see Figure [S]). 

The results here substantially strengthen the growing 
consensus that early galaxies are the key to reionization. 
These results on cosmic reionization from current state- 
of-the-art microwave background probes (PC15) com¬ 
bined with deep probes of faint z > 6 galaxies provided 
by the Hubble Space Telescope (e.g., Bouwens et al. 2011; 
McLure et al. 2013; R13; B15; Atek et al. 2015) have 
demonstrated the power of combining data from a wide 
range of major missions, and highlight the upcoming op¬ 
portunities with JWST for probing deeply into the reion¬ 
ization epoch. 
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APPENDIX 

A. PRODUCTION EFFICIENCY FOR LYMAN CONTINUUM PHOTONS PER UNIT UV LUMINOSITY 

In this section, we determine how the production efficiency per unit UV luminosity ^jon depends on the C/V-continuum 
slope P (where /? is defined such that f\ oc A^). Given that observations only allow for a constraint on the product of the 
two unknowns pionfesc (and perhaps more generally the product of the four unknowns Xionfescf{Miim){CHii/^)~^'^- 
see Eq. [5] from §4.2), choice of a physically plausible value for ^ion ensures that the fesc we infer is also of more 
relevance. 

Following a similar procedure to that executed in R13, we can estimate the efficiency parameter ^ion by considering 
a variety of different ages, metallicities, and dust content for star-forming galaxies at z > 6. For convenience, the star- 
formation rate of our model galaxies is assumed to be constant. The Chariot & Fall (2000) dust model is assumed, 
and we also leverage the spectral synthesis models of Bruzual & Chariot (2003). The CV-continuum slope P is derived 
from the model spectra over the wavelength 1700A and 2200A, consistent with the position of the broadbands used 
to derive P for z ^ 1 galaxies (e.g., Bouwens et al. 2010, 2014; Dunlop et al. 2013). The conversion factors and /3’s 
computed for many different model spectra are presented in Figure [TUI 

Fitting the envelope of derived conversion factors ^ion versus P, we find io25 i3-i.i(/3-i-2) s-i/(ergs s“^Hz“^) for 
P < —2 and io25 i3-o.6(/3-i-2) s-iy'(ergs s“^Hz“^) for P > —2, with an approximate width of this distribution of ±0.125 
in logio^on- 

B. EVOLUTION INFERRED FOR COSMIC IONIZING EMISSIVITY EXCLUDING CONSTRAINTS BASED ON THE 

PREVALENCE OF LYa EMISSION IN Z > 6 GALAXIES 

As discussed extensively in the text and elsewhere (e.g., Mesinger et al. 2015; Choudhury et al. 2015), large uncer¬ 
tainties exist in the use of Lya emission from galaxies to constrain the ionization state of the universe at z > 6. Not 
only do the observational inferences depend on the intrinsic velocity offset of Lya emission (e.g., Mesinger et al. 2015; 
Stark et al. 2015), but it is dependent on details of the radiative transfer, inhomogeneities in the IGM, the ionizing 
background, all of which rely on high-quality simulations of the reionization process (and where it can be challenging 
to include all of the relevant physics). 

As such, it is perhaps useful to derive the evolution in the cosmic ionizing emissivity without relying on observational 
constraints that involve the prevalence of Lya emission in z ^ 7 and z ~ 8 galaxies. We therefore repeat the exercise 
we performed in §4.1 but excluding constraints 3-4 from Table [T] 

Our derived constraints on the evolution of the ionizing emissivity (68% and 95% confidence intervals) are shown 
in the left panel of Figure [TTJ The equivalent (dlog,o yon/d-z)z=8j after marginalizing over parameter space and 
weighting by |Vt| (a flat prior in r: see §4.1), is —O.OO^o'is (Table [3]), which is 0.5cr shallower than in our fiducial 
analysis {dotted black line gives the 68% confidence intervals) where we fold in the constraints from the observed Lya 
fraction. The right panel in Figure [TT] compares these emissivity results with the evolution of the UV luminosity 
density (see also Figure [6|). 

Here again the evolution we infer for the ionizing emissivity is in excellent agreement with the evolution found in the 
UV luminosity density (B15). As in our fiducial analysis, the similarity of the evolution found for these two quantities 
suggests that star-forming galaxies provide the ionizing photons which drive the reionization of the universe. The 
exercise we present in this appendix shows that this conclusion does not depend on our use of constraints that involve 
the prevalence of Lya emission from z > 6 galaxies. 
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Fig. 10.— A determination of how the production efficiency ,Jion of Lyman-continuum photons per unit UV luminosity at I 6 OOA depends 
on the f/V-continuum slope /3. These efficiencies are calculated from the Bruzual & Chariot (2003) spectral synthesis library assuming a 
constant star formation rate. We also adopt three different metallicities (0.2 Zq, OAZq, and Zq) and a wide range in ages (0.1 Myr to 10 
Gyr). Both the case of no dust content and Ay = 0.1/0.2 (Chariot Sz Fall 2000) is considered, as indicated on this figure. ^ is computed 
over the spectral range 1700A to 2200A. The shaded envelopes indicate the approximate dependence of ^lon on /3- For completeness the 
full range of dusty and dust-free tracks are shown, even up to solar metallicity. We recognize that solar metallicities may be unlikely at 
z > 6 where results (e.g., Bouwens et al. 2014; Dunlop et al. 2013) suggest that the dust content is very low. 

C. COMPARISON BETWEEN THE INTEGRATED r’S AND REIONIZATION HISTORIES IMPLIED BY PLANCK AND THE 

WMAP 9-YEAR RESULTS 

For context and to illustrate the gains from Planck, we also Figure [T^ showing the filling factor evolution of ionized 
hydrogen, Qhii, and integrated Thomson optical depths r’s as inferred from the WMAP 9-year r measurement and 
assuming reionization is complete at z = 5.9-6.5 {magenta hatched region: 68% conhdence). 
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Fig. 11.— (left) 68% and 95% confidence intervals (red-shaded region) on the evolution of the cosmic ionizing emissivity (assuming 
Chii = 3) supposing that reionization is complete aX z = 5.9-6.5, utilizing the latest Planck optical depth measurements (PL15), and 
enforcing continuity with inferred ionizing emissivity aX z = 4.75 (constraints 1, 2, and 5 from Tabled and excluding constraints based 
on the prevalence of Lya emission in galaxies (Appendix B). The dotted black lines demarcate the 68% confidence regions for our fiducial 
determination. See the caption to Figure [6] for a description of other symbols and regions, (right) Comparison of the evolution of the 
inferred ionizing emissivity with the evolution of the UV luminosity density (black-hatched region) adopting some redshift-independent 
conversion factor log^Q fesdion = 24.53. As before, the evolution of the ionizing emissivity is very similar to that found for the luminosity 
density of galaxies in the rest-frame UV. While this again suggests that star-forming galaxies drive the reionization of the universe, this 
figure demonstrates that this conclusion does not depend on whether or not we make use of constraints that involve Lyo; emission from 
z > 6 galaxies. 
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Fig. 12.— Comparison of the constraints on the Thomson optical depth r provided by the Planck 3-year results (PC15: cross- 

hatched black region) and the cumulative optical depths preferred at 68% and 95% confidence by our models of the ionizing emissivity 
evolution including new constraints from Planck. The magenta-hatched region show the constraints that would be provided by the WMAP 
results (68% confidence), {right) Shown are constraints on the filling factor of ionized hydrogen Qhii a function of redshift and the 
evolution preferred at 68% and 95% confidence by our models of the ionizing emissivity evolution. The magenta-hatched region indicates 
the range of Qhii allowed at 68% confidence for the WMAP 9-year r measurement (0.089 ib 0.014). Similar to Figure!^ 










